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The rotational temperature of polar molecular ions in Coulomb crystals 
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With MgH + ions as a test case, we investigate to what extent the rotational motion of smaller 
polar molecular ions sympathetically cooled into Coulomb crystals in linear Paul traps couples to the 
translational motions of the ion ensemble. By comparing results obtained from rotational resonance- 
enhanced multiphoton photo-dissociation experiments with data from theoretical simulations, we 
conclude that the effective rotational temperature exceeds the translational temperature (< 100 
mK) by more than two orders of magnitude, indicating a very weak coupling. 
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Molecules which are translationally cold, spatially lo- 
calized and internally prepared in specific states consti- 
tute ideal targets for a large variety of investigations, in- 
cluding high resolution spectroscopy [l|, coherently con- 
trolled chemistry [2| , state-specific reaction studies 
and molecular internal state control experiments 

For atoms such targets have for years been available 
by combining laser cooling techniques with trapping and, 
perhaps most spectacularly, been exploited in quantum 
logics experiments 0,0- It is, however, extremely diffi- 
cult to adopt the same techniques for molecules due to 
the general lack of closed optical transitions. For polar 
molecules a further complication arises since the rovibra- 
tional degrees of freedom couple rather efficient to the 
black body radiation (BBR) of the surroundings. Re- 
cently, various methods that at least partially approach 
the ideal situation have, however, been demonstrated for 
neutral molecules. These include photo-association of 
cold atoms 13, buffer gas cooling of molecules in mag- 
netic traps [9j , and control of the motion of molecules by 
either by static 0] or optical electrical fields [ill ]. 

Internal cold, but not spatial localized molecular ions 
have been produced either by photo-process es 1121 ll.1l . 
by buffer-gas cooling in RF multi-pole traps |l4 . Il5| or 
by state selective recombination llfjj. Recently, by sym- 
pathetical cooling to translational temperatures in the 
mK-range in RF ion traps, spatially localized molecular 
ions in Coulomb crystals have been reported [HI flif . 

In this Letter, we investigate whether there is a suf- 
ficient coupling between the internal degrees of freedom 
of polar molecular ions and the external motion of the 
ions in a Coulomb crystal to achieve a significant sym- 
pathetic cooling of the internal states of the molecular 
ions. The strategy is to apply Resonance-Enhanced Mul- 
tiphoton Photo-Dissociation (REMPD) 0] as a tool to 
obtain knowledge of the populations in the various rota- 
tional states and hence determine the effective rotational 
temperature. One could most likely expect a measurable 
cooling of the internal degrees of freedom of the molecular 
ions due coupling of the dipole moment of the molecules 
to the motion of the charges of the ions in the Coulomb 
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FIG. 1: The four lowest spin singlet potential curves for 
24 MgH + including various symbols and processes discussed 
in the text. 



crystal. However, since the typical frequencies of the nor- 
mal modes of the Coulomb crystal will be ~ 10 6 Hz while 
the rovibrational frequencies will be ~ 10 11 Hz - ~ 10 13 
Hz, no strong resonant coupling is expected Ref. |2(i|. 

An ensemble of 24 MgH + molecules, which is trans- 
lational sympathetically cooled by laser cooled 24 Mg + 
atomic ions in a linear Paul trap has been chosen as a 
test case for many reasons. The 24 MgH + ions are easy to 
produce and sympathetically cool |l7| , at room tempera- 
ture 99% will be in the vibrational ground state, precise 
spectroscopic data are available for 24 MgH + [il|, and fi- 
nally in order to compare experimental results with the- 
oretical simulations, the potential curves and permanent 
as well as transition dipole moments of 24 MgH + can be 
calculated with good precision [23, |23j . 

In Fig. ^ the energy potential curves of the four lowest 
lying spin singlet electronic states of 24 MgH + are shown 
together with an indication of states and processes of 
relevance to the investigations. The resonant part of 
the two-photon REMPD process consists of a rotational 
state-selective excitation from a state \v = 0, J)x in the 
X 1 E potential to a state \v" = 0, J + 1}a in the A 1 !] po- 



tential by a laser field e. The letters in the ket-notation 
represent the vibrational and rotational quantum num- 
bers, respectively. From the \v" = 0, J + 1)a state, the 
molecule can dissociate nearly exclusively into Mg and 
H + via the C potential curve by absorbing another pho- 
ton from the e field. However, with the field strengths 
applied in our experiments, the excited molecule primar- 
ily spontaneously decays back to a vibrational excited 
state in the electronic ground state potential X 1 !! with 
a rate Te ~ 3 x 10 8 s _1 . Subsequently, further sponta- 
neous emissions, stimulated processes due to the coupling 
to the BBR field and transitions induced by an effective 
coupling to the external motion of the ions take place. 
In Fig. I, these various mechanisms are in brief denoted 
by r„j. Eventually, the molecule will reach the state 
\v = 0, J)x from where it again can couple strongly to 
the laser field e. 

The experimental setup has previously been described 
in Refs. |2J,|25j. In short, 24 Mg + ions are produced by 
off-resonance photo-ionization by the pulsed laser field e 
indicated in Fig. 1 and trapped in a linear Paul RF trap 
[2^ | situated in a UHV chamber at background pressure 
of ~ 1 x 10~ 10 torr. The produced 24 Mg + ions are laser- 
cooled on the 3s 2 Si/2-3p 2 P3/2 transition by two cw laser 
beams counter-propagating along the center-axis of the 
trap to temperatures below 10 mK where they form a 
Coulomb crystal [24J. The two-dimensional projection 
images of such crystals are obtained by monitoring the 
fluorescence by a CCD camera. A typical projection im- 
age of a pure 24 Mg + crystal is shown in Fig. Efa). The 
shape of the real crystal is spheroidal with rotational 
symmetry around the dotted line. 

The 24 MgH + molecular ions are formed through pho- 
tochemical reactions between 24 Mg + ions excited to the 
3p 2 P 3 / 2 state and H2 molecules leaked in the chamber 
at partial b ackg round pressure of 5 x 1CP 10 torr for a 
few minutes [l7|. The produced 24 MgH + molecular ions 
arc within tens of milliseconds translational sympathet- 
ically cooled by the remaining 24 Mg + to a temperature 
below 100 mK, and a two-component Coulomb crystal is 
formed 01- As seen m Fig. IIP 3 ), due to dynamical ra- 
dial confinement of the ions, the non-fluorescing 24 MgH + 
molecular ions and the coolant 24 Mg + ions segregate ra- 
dially with the lighter mass ions, the 24 Mg + ions, closest 
to the rotational symmetry axis [24j . The amount of each 
ion species in the crystal can be deduced from the fact 
that each species in the trap will have a uniform density 
inversely proportional to its mass and that for crystals 
containing two species with a small difference in mass, 
the outer boundary will have a near spheroidal shape 
[25H26I |2?| . By determining the volume of the fluorescing 
region in Fig. ffib). we can conclude that this particular 
crystal has a relative 24 MgH + content of ~ 65%. 

The laser field e presented in Fig. ^ originates from 
a frequency doubled dye laser pumped by a Nd:YAG 
laser. The pulse repetition rate is 20 Hz, and the dura- 



2 




FIG. 2: (a) 2D projection-image of a pure Coulomb crystal of 
approximately 3000 fluorescing 24 Mg + ions. The length of the 
crystal is 1070 /im. (b) A two-component crystal, consisting 
of ~ 35% 24 Mg+ and ~ 65% 24 MgH+ ions (non-fluorescing), 
obtained by letting 24 Mg + ions in the crystal in (a) react with 
a gas of H2. (c) A nearly pure 24 Mg + crystal obtained after 
ns-laser pulses have been applied for 60 s to dissociate the 
24 MgH + molecules in the crystal shown in (b). The dashed 
ellipse indicates the position of the outer surface of the initial 
crystal in (a). The dotted line in (a) indicates the axis of 
rotational symmetry. 



tion of each pulse is approximately 5 ns. The waist of the 
field at the trap center is 500 /im, and the peak intensity 
was typically of the order of I pea k ~ 10 7 W/cm 2 , corre- 
sponding to a peak power-broadening of the \v = 0, J)x _ 
\v" = 0, J + 1)a transitions of ~ 0.3 cm -1 . The spectral 
width of each pulse was ~ 0.25 cm -1 , and the center 
frequency was calibrated to better than 0.1 cm -1 . Dur- 
ing the experiments, the center frequency was scanned 
0.95 cm^ 1 at 1 Hz in order to obtain results which are 
independent of slow frequency drifts in the laser sys- 
tem. The resonance frequencies of the X 1 S(i/ = 0, J)- 
A 1 S(^" = 0, J + 1) transitions, which are used in our 
REMPD scheme, are known with a precision of 0.05 cm -1 
fromRef. Hj. 

Before the REMPD lasers are applied, the produced 
24 MgH + ions are allowed to internally equilibrate for 
about one minute, a time long enough to ensure that the 
produced molecular ions arc in the vibrational ground 
state and that any effective cooling of the rotational de- 
grees of freedom will be revealed in the experiments. In 
the experiments, crystals like the one shown in Fig. 2(b) 
are exposed to the laser pulses while CCD images are be- 
ing recorded at a rate of 6 Hz. After a typical exposure 
time of tens of seconds, all molecules have been dissoci- 
ated, and a nearly pure crystal of 24 Mg + is left in the 
trap as shown in Fig. Hfc). Before a new REMPD ex- 
periment is carried out, the trap is reloaded with 24 Mg + 
ions and through reactions, Coulomb crystals with nearly 
the same number of 24 Mg + and 24 MgH + ions are created 
each time. 

In Fig. |31 the normalized amount of 24 MgH + molecular 
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FIG. 3: The amount of 24 MgH+ ions relative to the initial 
content of 24 MgH + in the crystal as a function of the time 
in which the 24 MgH + ions are exposured to the dissociation 
laser pulses. The triangles and circles represent data for the 
R(l) and R(6) resonances, respectively (see text). In this case 
Wk = 1.9 x 10 7 W/cm 2 . 




Cooling rate r c [s 

FIG. 4: The steady-state rotational temperature of 24 MgH + 
as a function of the effective cooling rate F c (see text for 
detailed definitions). The insert shows the population dis- 
tributions for three different rotational temperatures: White: 
T rot =300 K, Grey: T rot =120 K, and Black: T rot =7 K. The 
colors of the bars correlates to the color of the three squares 
in the main plot. 



ions are shown for two out of about ten REMPD experi- 
ments where the e field has been tuned to the resonance 
of the R-branch transitions R(l): \v = 0, J = l)x~ \v" = 
0, J = 2} A and R(6): \v = 0, J = 6) X "K = 0,J = 7) A , 
respectively. 

In order to derive an estimate for the rotational tem- 
perature from the outcome of the REMPD experiments, 
simulations including various strengths of a effective sym- 
pathetic rotational cooling have been performed. In the 
model, this cooling has been accounted for by introduc- 
ing an artificial decay rate from the rotational states 
J to J-l proportional to J. This corresponds roughly 
to have an effective cooling rate T c defined through 
dE rot /dt = —T c E rot , where E rot is the average rota- 
tional energy. In Fig. 4, the simulated steady-state rota- 
tional temperature of the 24 MgH + ions, defined as T rot 
= E rot /ks , is presented as a function of the cooling rate 
T c . In the absence of sympathetic cooling, the rotational 
distribution is observed to equilibrate to the temperature 
of the BBR field in the trap region (~ 300 K). It is, how- 
ever, clear that only a very modest cooling rate r c (~ 0.1 
s" 1 ) will change this situation dramatically (see Fig. 4). 
The steady-state rotational distributions are the stating 
points for all simulations. 

When simulating the time evolution of the photo- 
dissociation process, the time scale is broken into two 
types of intervals defined by the presence or absence 
of the laser pulses. During a laser pulse, the effect of 
laser driven transitions between the various rovibrational 
states of the X X E and A X E potential curves is taken into 
account by at any time during the pulse assuming that 
the various populations are in steady-state with the in- 
stantaneous laser field strength and detunings. Based 



on the instantaneous populations in the A 1 E potential 
curve, the dissociation rate via the C^E potential curve 
is determined. During the laser pulse, spontaneous emis- 
sion from the populated A 1 !] states to the various rovi- 
brational states of the X 1 E potential curve is included 
as well, while rovibrational transitions within the poten- 
tial curves are neglected due to the short time scale of 
the pulses. Between consecutive pulses population re- 
distribution processes among the rovibrational states in 
the electronic ground potential X 1 E are simulated taking 
into account a effective sympathetic cooling, spontaneous 
emissions and interactions with the BBR field. Due to 
the complexity of the simulations, only the main ideas 
have been considered here, but the details will be pre- 
sented elsewhere Ref. |2fjj . 

The ratios between the simulated dissociation rates for 
REMPD processes via the R(l) and R(6) resonances are 
presented in Fig. |SJas a function of the effective rotational 
temperature obtained by including a certain cooling rate 
r c in the calculations. The dissociation rates have been 
extracted 10 s after initiation of the photo-dissociation 
process where the various populations and the dissocia- 
tion rate are simulated to be in steady state [2(i|. The 
chosen dissociation ratio provides a simple experimental 
way of distinguishing between low (< 1 K) and high (> 
100 K) rotational temperatures T rot as the ratio changes 
drastically for temperatures below ~ 50 K. 

The ratio of dissociation rates is fortunately nearly in- 
sensitive to intensity variations. Even a factor of four in 
the actual peak intensity does not significantly change 
the ratio. This indicates that a smaller error in cali- 
bration of the power meter used in the experiments or 
smaller uncertainties in the laser beam positioning should 
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FIG. 5: The ratios between the simulated dissociation rates 
for REMPD processes via the R(l) and R(6) resonances. The 
values of T rot are directly connected to a cooling rate F c in- 
cluded in the simulations (presented in Fig. 4). The curves 
represent three different peak-intensities of the laser pulses, 
whereas the hatched area represents a conservative estimated 
range derived from experimental data as those presented in 
Fig. 3. 



not affect the measured ratio. A conservative estimated 
interval of this ratio that is in agreement with a series 
of experimental data as those presented in Fig. 3 is in- 
dicated in Fig. 5 as the hatched area. From Fig. 5 it 
is hence reasonable to conclude that the rotational tem- 
perature of the 24 MgH + ions is higher than 120 K, and 
likely as high as room temperature, which is the same ro- 
tational temperature found in storage ring experiments 
with very weak coupling between the individual ions |28l | . 
In terms of the rotational cooling rate r c , we can corre- 
spondingly conclude that it must be smaller than ^0.1 
s — 1 . Since the permanent dipole moment of 24 MgH + (3.6 
Debye) is not particularly small, and its rotational fre- 
quencies are not extraordinary, we expect that molecular 
ions in Coulomb crystals do not generally become rota- 
tionally cold. Also since the typical trapping potentials of 
Penning traps are typically similar to those of RF traps, 
no gain with respect to rotational cooling in such traps 
is expected. 

When experimenting with polar molecular ions at 
room temperature, it hence seems necessary to imple- 
ment a kind of active rotational cooling in order to 
achieve state specific molecular ion targets [13, 13, ■ 
Alternatively, one might choose to cool the trap envi- 
ronment to cryogenic temperatures |l4l llij . For non- 
polar molecular ions, the coupling between the internal 
degrees of freedom and the modes of the Coulomb crys- 
tal as well as the BBR field is smaller than for polar 
molecules. Hence, such ions produced by, e.g., a state- 
selective REMPI-process would possibly be able to reside 
in this state through sympathetic cooling, and could lead 
to spatial localized state-specific targets. 



In conclusion, with 24 MgH + ions as the test case, we 
have proven that the rotational motion of smaller polar 
molecular ions sympathetically cooled into a Coulomb 
crystal does not become significantly cooled, but very 
likely stays in nearly thermal equilibrium with the sur- 
rounding environment. 
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